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Quantum imaging promises increased imaging performance over classical protocols. However, there are a num- 
ber of aspects of quantum imaging that are not well understood. In particular, it has so far been unknown how 
to compare classical and quantum imaging procedures. Here, we consider classical and quantum imaging in a 
single theoretical framework and present general fundamental limits on the resolution and the deposition rate 
for classical and quantum imaging. The resolution can be estimated from the image itself. We present a utility 
function that allows us to compare imaging protocols in a wide range of applications. 
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Imaging is an important technological tool in many disci- 
plines including astronomy, biomedical research, nanotech- 
nology and basic research. The wave nature of light limits the 
resolution and contrast that can be achieved in classical imag- 
ing techniques. On the other hand, quantum entanglement may 
offer some improvement over these classical limits, leading to 
the subject of quantum imaging. The best known quantum 
imaging protocols are two-photon microscopy [1-3] and spec- 
troscopy [4], quantum holography [5], quantum lithography 
[6], and quantum illumination [7]. The subject of quantum 
imaging is closely related to quantum metrology [8], which is 
quite well understood [9, 10]. However, there are some impor- 
tant differences. If we restrict our discussion to quantum op- 
tics, quantum metrology is first and foremost concerned with 
photon statistics, which is governed almost entirely by transfor- 
mations of the mode operators of the quantum field. Quantum 
imaging, on the other hand, is intimately related to both the 
operator transformations and the specific form of the (classi- 
cal) mode functions that encode the spatial properties of the 
field. This has made it difficult to identify the essential quan- 
tum mechanical behavior in imaging protocols. 

There are two reasons why it is important to understand 
the precise distinction between classical and quantum imag- 
ing. First, it can help identify new methods for improved imag- 
ing, potentially leading to new technologies. Second, it may 
reveal a fundamental aspect of physics that has hitherto re- 
mained elusive, namely what (if any) makes quantum optics 
more powerful for imaging than classical optics. The obvious 
answer to this question — quantum entanglement — has already 
been proved false to some extent [11]. While entanglement 
is likely necessary for an improved imaging procedure, it is 
certainly not sufficient. This is reminiscent of quantum com- 
puting, where it was shown that entanglement is necessary but 
not sufficient for obtaining the promised exponential speed-up 
over classical computing [12]. 

Here we develop a theoretical framework that allows us 
to study classical and quantum imaging in a unified fashion, 
which carries over tools and techniques from the mathematical 
theory of metrology to prove fundamental limits. This is fun- 
damentally and conceptually different from studying parame- 
ter estimation systems that happen to use light as the probing 




FIG. 1. General imaging procedure. The state of light p(0) contains 
information about the object O and the imaging system L. The sub- 
strate S acts as a measurement of p(6) and is generally described by 
a povm. This model includes multi-photon quantum imaging proce- 
dures. 



mechanism (see e.g., [13]). Our methodology is based on the 
statistical distance between two probability distributions that 
characterise two images. We find that the distinguishability 
of the probability distributions provides a natural definition of 
the imaging resolution in terms of a Cramer-Rao bound. We 
apply this bound to several imaging procedures and show that 
it gives the desired result. We define a utility function that 
provides a general metric for imaging procedures in practical 
applications. 

The most general imaging procedure is shown in Fig. 1. 
Light emitted or reflected by an object O is optically processed 
by an imaging system L and recorded by a substrate S. Typ- 
ically, the substrate has a granular spatial structure so we can 
break it up into N systems indicated by their position x. If 
the width of such a "pixel" is £, then the total length of the 
substrate is L = NL Moreover, we assume that the inten- 
sity is granular (which is consistent with the use of photons at 
the fundamental level), and the values of the image intensity 
at a pixel are given by i^. The average recorded intensity at 
pixel x is then given by I(x) = J2k ^kPk(x), where Pk(x) 
is the probability of measuring at pixel x. In quantum me- 
chanics, this is given by the Born rule Pk(x) = Tr[pEk(x)], 
with p the quantum state of light incident on the substrate and 
Ek(x) the measurement operator for outcome ik at x. The set 
of Ek (x) form a positive operator- valued measure (povm) for 
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each pixel x. The image I(x) may be related in a simple man- 
ner to the intensity of the light that hits the substrate, or it may 
have a more complicated relationship involving multi-photon 
absorption, etc. For monochromatic images, I(x) will be a 
single-valued function indicating the brightness of the image 
at x. Note that the pixels do not need to be placed side-by- side 
(see, e.g. Ref. [14]). For simplicity we will limit our discus- 
sion to one-dimensional substrates in the remainder; however 
all our results generalise easily. 

We can construct the image observable I(x) that is mea- 
sured by the substrate. We define I(x) = (I(x)) such that 
H x ) — J2k T^[pEk(x)] = Tr [pi (x)], and we therefore 
find that I(x) = J2 k ik Ek{x). In the simplest case where the 
substrate records the number of photons k, the povm elements 
are Ek{x) = \k) x (k\ and ik = k. More complicated povms 
can incorporate for example photon losses, a saturation point, 
bleeding from adjacent points, or multi-photon interference. 
In typical imaging protocols the measurement at each pixel is 
diagonal in the Fock basis, and the properties of the substrate 
are characterised by the povm 



Ek(x) = ^2qk(x\n)\n) (n\ 
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where qk(x\n) is the conditional probability of finding mea- 
surement outcome k at pixel x given the photon number distri- 
bution n = (774, . . . , tin) over A 7 " pixels, andj^ qk(x\n) = 1. 
We identify the recorded image I(x) with an unnormalized 
probability distribution. After normalisation we define 
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where Iq = I(x) is the image normalisation. 

In imaging, we are often interested in distinguishing two 
possibilities, which we can parametrize by a continuous num- 
ber 6, for example the position of a light source. Hence, 
the probability distribution Pr(x\6), and therefore the image, 
will depend on 0. We now have two probability distributions, 
Pt(x\6 = 0) and Pr(x\0). Distinguishing two images is then 
possible if we can distinguish the two probability distributions. 
One method to achieve this is to define the infinitesimal statis- 
tical distance [15] 
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[Pr(x\d0) -Pr(>|0)] 2 
Pr(>|0) 
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which can be formally integrated to a function s(0). We note, 
however, that in principle any method for distinguishing prob- 
ability distributions can be used. Using a Taylor expansion of 
s (0) up to first order in 50, we can write 



5s = s(S0) - 5(0) 



ds 
dO 



(4) 
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The two probability distributions are distinguishable when 
the statistical distance between them exceeds a constant value 



(Ss) 2 > 1, which we can rewrite as 



(5s) 2 = F(0)(50) 2 > 1 with F(0) 



• (5) 



The quantity F(6) is known as the Fisher information [16] 
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and is a measure of the information we can extract about 6 
given the probability distribution Vx{x\9). When we use the 
form of Pr(a:|0) in Eq. (2) we find 
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where 3q = d/d6. If 7rj is independent of 6 the last term van- 
ishes. The Fisher information is therefore calculated directly 
from the recorded image, without the need for modelling the 
povm {Ek(x)} of a complicated substrate. 

We can translate the distinguishibility criterion in Eq. (5) 
using Fq = F(0) into a resolution limit for the parameter 
and find that 
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This is a single- shot Cramer-Rao bound for the imaging reso- 
lution of 0, which is achievable in the limit of a well-exposed 
image. This result is readily extended to multiple parameters 
= . . . , 0d) using well-known techniques [16]. 

There is a key difference between the theory of imaging pre- 
sented here and the theory of quantum metrology [17]. In the 
latter, the Cramer-Rao inequality gives a bound on the mean 
square error in the parameter as a function of the number 
of independently repeated measurements [13]. By contrast, 
in quantum imaging the Cramer-Rao bound is based on the 
full probability distribution Py(x\0), and hence assumes that 
the image is sufficiently well-exposed such that Eq. (2) holds. 
Faint images can be modelled phenomenologically by adding 
a (uniform) probability distribution to Eq. (2). 

Next, we demonstrate that this technique gives the correct 
results for some known imaging procedures. First, we con- 
sider the double slit experiment, in which the image is the fa- 
miliar interference pattern in the far field [18]. We parame- 
terise the slit separation by and have the numerical aperture 
A capture a fraction of the interference pattern. To calculate 
the Fisher information we use the full form of Eq. (7) since 
Io will be heavily dependent on 0. We then identify that the 
resolution criterion is that the uncertainty 50 is smaller than 
the estimate 0. From this we derive that two slits are resolved 
when 2 F{0) > 1. We then calculate the resolution limit to be 
#min > 0.369A/A When we compare this to the Abbe limit 
#min > 0.5 A/ A we find a slightly better resolution. This is 
because the Abbe limit is a conservative estimate that does not 
take into account the extra information about in the slopes of 
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substrate 



FIG. 2. Optical lithography. Two nearly counter-propagating beams 
a and b produce an interference pattern on the substrate. The phase 
shift causes a shift in the image. 



the intensity pattern. The Fisher information is very sensitive 
to these slopes. When we calculate the resolution for the third- 
order photon interference slit imaging experiment by Oppel et 
al., [14], we find the factor 2 improvement, as expected. 

Second, we consider classical and quantum photo-litho- 
graphy [6] shown in Fig. 2. In the classical regime the in- 
put state can be modelled by the single photon superposition 
state = (|l,0) a6 + e i9 \0,l) ah )/V2, where = k£Ax 
shifts the pattern along the substrate (k = 2tt/X is the wave 
vector). We give the povm of a general M-photon absorbtion 
substrate including a detection efficiency rj. The conditional 
probabilities in Eq. (1) are qo(x\n x ) = 1 for n x < M and 
Qo(x\n x ) = (l-r]) nx forn x > M\ qi(x\n x ) = l-q (x\n x ). 
The Fisher information for M = 1 is then calculated as 

F w = ^E sin2 ( Kte+ 9 w 

If the substrate is much larger than the wavelength of the light, 
then Fo ~ 1 (see Fig. 3a) and the Cramer-Rao bound is SO > 1, 
or £Sx > A/27T. This imaging resolution is slightly better than 
the Rayleigh criterion £Sx > A/4 due to the slope sensitivity 
of the Fisher information. 

For M-photon quantum lithography the Fisher information 
is given by 

When kIN ^> 1 the Fisher information is Fo ~ M 2 , yielding 
SO > M _1 . The imaging resolution iSx > A/27rMis slightly 
better than iSx > A/4M found by Boto et al [6]. In both the 
classical and quantum case the resolution is independent of the 
efficiency r]. 

Finally, we consider imaging the position of a dot of 
coherent light with a Gaussian intensity profile on a one- 
dimensional substrate. We include a saturation limit of the 
substrate and pixel bleeding due to interactions between adja- 
cent pixels. For simplicity the state of light on the substrate 
is given by \a) = (£) x=1 \ot x )x with \a x ) x a coherent state of 
amplitude a x = ao exp[— (x — xq) 2 /2a 2 ] at pixel x. There is 
no coherence between adjacent pixels. Assuming a <C L and 
xq towards the centre of the substrate (so DqIq ~ 0), the Fisher 
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FIG. 3. (color online) Qualitative behaviour of the Fisher informa- 
tion: (a) Fo in Eq. (9) as a function of the number of pixels N tends 
to 1 for large N; (b) Fo for the displacement of a Gaussian dot with 
width a (in units of pixel size £). The vertical dotted line indicates 
where a = £, the size of one pixel; (c) F for a Gaussian dot as a 
function of the peak amplitude ao (in arbitrary units). The vertical 
doted line indicates the saturation point at ao = 5, and a — 70£. The 
saturation point has only a small effect, since Fo is sensitive mostly 
to intensity slopes; (d) The ratio of the Fisher information with and 
without bleeding Fb /Fo as a function of the average bleeding distance 
b (with a 40£ and a 70). 

information for perfect photodetection at each pixel becomes 

F(fl)=4 j]b£^. (id 

^ Iq cr z 

X 

The Fisher information is shown as a function of the standard 
deviation a in Fig. 3b. When a falls below one pixel width 
the Fisher information drops rapidly to zero. This is expected 
since below one pixel width we are no longer able to detect a 
correspondingly small shift in 0. Above one pixel width the 
graph follows the form of F = 2/<j 2 , which gives the resolu- 
tion SO > a/V2. 

Taking into account the saturation limit S leads to condi- 
tional probabilities in the povm given by 

qk{Xln)= \ET^n M if * = 5. (12) 

For the displacement of a Gaussian dot this gives the Fisher 
information shown in Fig. 3c. Due to a uniform noise floor Fq 
rises sharply for small ao and falls away at the saturation point. 
The noise is numerically necessary to avoid near- singularities. 
The saturation truncates relevant statistical information, but 
the main contribution to Fo comes from the steep ascent re- 
gions of the Gaussian, rather than the crown. Therefore, the 
overall impact on the Fisher information is relatively small for 
this particular imaging process. In other situations, the satura- 
tion limit can have severe implications for the resolution. 

When there exists an interaction between adjacent pixels, 
a signal at pixel position x is recorded by its neighbours 
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with some probability. This is called bleeding, and the con- 
ditional probabilities of Eq. (1) are given by qk(x\n) = 
Yl x ' ^fc,n cc /Pi > ((i), where Pr(d) is the probability of a signal 
bleeding a distance d = \x' — x\ . For a Poisson distribution 
the average distance is characterised by the bleeding parameter 
b. The qualitative effect on the Fisher information is shown in 
Fig. 3d. 

Given the Cramer-Rao bound for the imaging resolution, 
we would like to find a physical limit to the Fisher informa- 
tion that in turn would give a fundamental limit to all imaging 
procedures. The Fisher information is bounded by F(0) < 
4(AK) 2 /h 2 [19], where AK is the variance of K, the gen- 
erator of translations in 6 such that p(8) = Up(0)U^ with 
U = ex.p(—i6K/h). However, this relation does not include 
any information about the substrate, since it is a general bound 
on any possible metrological system including ideal imaging 
procedures. Another important question therefore concerns 
the best possible resolution over all imaging systems and sub- 
strates given a particular constraint on the substrate. For ex- 
ample, what is the best imaging resolution given a certain sat- 
uration limit S, or a given amount of pixel bleeding bl This is 
currently an open question. 

As pointed out by Tsang [20] and Steuernagel [21], the prac- 
ticality of some quantum imaging procedures such as quan- 
tum lithography [6] is severely affected by the efficiency of the 
measurement procedure. In the case of M-photon lithography 
the rate of the M-photon detection events at pixel x yields our 
image I(x). This is calculated as 
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The constant of proportionality disappears in the normalisa- 
tion procedure of Eq. (2), and hence does not contribute to 
the Fisher information. However, it is clear that the deposition 
rate of M-photon detection events falls off exponentially. Any 
practical figure of merit must therefore take into account both 
the resolution and the deposition rate of the imaging proce- 
dure. We define a utility function U for an imaging procedure 
as U = FqD c , which should be maximised. Here, < D < 1 
is the deposition rate and c > is a cost function that can 
take into account the resource cost of each attempted detection 
event. This may reflect both physical and economical aspects 
of the imaging procedure. For example, in the manufacturing 
of microchips speed is an important factor and the cost func- 
tion should be suitably weighted to reflect this. The deposition 
rate can be defined based on the povm element Eq that indi- 
cates a failure to detect anything. If Eq (x) is the povm element 
for not detecting anything at pixel x, then E = ® x=1 Eq(x) 
and the deposition rate is given by 



N 



X = l 



(14) 



For M-photon lithography, D oc A/'-( M - 1 ). For M > 1 this 



becomes a function that decreases polynomially with the num- 
ber of pixels N. 

In conclusion, we have presented a general method for es- 
tablishing the quality of classical and quantum imaging proce- 
dures based on the Fisher information. The resulting Cramer- 
Rao bounds give better estimates of the resolution than the 
traditional Abbe and Rayleigh criteria, and can also take into 
account the imperfections and edge effects of the substrate. 
Moreover, this technique can be used in imaging via higher- 
order photon correlations. The Fisher information can be cal- 
culated using a povm that models the substrate, or it can be 
inferred directly from the experimental data itself, circumvent- 
ing the need for complex mathematical modelling. We found 
that a saturation limit of the substrate has relatively small ef- 
fect on the imaging resolution. A utility function was defined 
that takes into account the resolution as well as the efficiency 
of creating the image. 
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